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1. Introduction 



An important part of our experimental understanding of hadron substructure is that hadrons in 
a medium are different to those in free space. The "EMC effect", the modification of the proton 
p2{x,Q^) structure function inside a nucleus (as first observed in 1983 by the European Muon 
Collaboration [[T]]) is one of the most famous examples of medium sensitivity. This effect has 
been studied in many hadronic models with varying degrees of success (see Ref. for a recent 
review) but an understanding from first principles is lacking. While the EMC effect and other such 
modifications of hadron properties are not a posteriori unexpected from the viewpoint of Quantum 
Chromodynamics (QCD), the complexity of QCD calculations involving nuclei has prevented the 
direct investigation of such effects. In QCD, we expect that medium modification is a ubiquitous 
feature of complex hadronic systems and, in these proceedings, we report on the study of a close 
analogue of the EMC effect, namely the modification of pion structure in the presence of a Bose 
condensed medium of pions. 

Since hadronic structure is a low energy consequence of QCD, the only tool with which to 
perform ab initio studies is lattice QCD (LQCD). Lattice QCD is formulated in Euclidean space, 
so physics defined on the light-cone, such as that embodied in the parton distributions and structure 
functions of deep inelastic scattering, is very difficult to address directly (an alternate approach is 
suggested in Ref. ^). However, through the Wilsonian operator product expansion (OPE), the 
Mellin moments, (x")/,(/i) = j\dxx"qh{x;iJ.), of the unpolarised quark distribution, ^/,(x;/i), in 
a hadron h correspond to the forward matrix elements of local operators {jJ, is the renormalisation 
scale). Here, we focus on the leading twist, unpolarised operators and have 

{h; p\4^-^"^\h; p) = {^)hP^^...p^^, (1.1) 
^{/ifl...M„}(^) = q{x)y^^^D^'' . ..D^"\{x) . (1.2) 

where {...} indicates symmetrisation of enclosed indices and subtraction of traces. The depen- 
dence of the various quantities on the renormalisation scale is suppressed for concision and q 
represents a particular flavour of quark field. The hadron h can be a proton or pion or it can be 
a more complex object such as a nucleus or a collection of mesons. In these proceedings, matrix 
elements of the n = \ operator are investigated in systems of up to twelve pions. 



2. Lattice Methods 



These matrix elements can be computed using the lattice approach by measuring two- and 
three- point correlation functions in the appropriate hadronic states. For clarity, we will consider 
the case of the up quark distribution of the 71+ {ud) in a medium of 71+ 's. Two point functions from 
a source location xq = (xo,fo), 



C„,(f,p) = (0 



nie'P'-'';r+(x, 



i=l X 



(^o)]' 



(2.1) 



where 71+ = uy^d, allow the energies of systems of m-pions to be determined from the dependence 
on Euclidean time. The total momentum of the m-pion state is p = YT=\Vi selected by the 
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summations over spatial sink locations (the individual p, are not quantum numbers). In the current 
context, we shall only consider p = 0. 

For fo = 0, the spectral decomposition of the correlator has the form 

where T is the temporal extent of the lattice and the ellipsis denotes excited states that are expo- 
nentially suppressed as t increases. The factors ~ |(£|(7r+)'"|m — represent the overlap of 
the m-pion interpolating operator onto {m — £) tt+'s going forward in time and i 71^ 's going back- 
ward in time (in the ground state, the momentum of the forward and backward going collections 

(m—£) (£) 

of pions separately vanish). Note that Z,„ = Z,n . In the limit of a large temporal extent of the 
lattice (r — )• oo), only the term with i = contributes, but at finite T (corresponding to non-zero 
temperature), thermal states in which some number of pions travels around the temporal boundary 
are important as we shall see below. 

Corresponding three point correlation functions allow the matrix elements in Eq. ( |1.1[ ) to be 
determined for the operator ^^"^ = 



Ct\r,t,p) = (0 



(=1 X 



. (2.3) 



Here, the operator is inserted at time-slice T injecting zero momentum. The spectral decomposition 
of the three point correlator for p = is 

= tjy"i)^n {i?^:\)e-^^^^^^^ (2.4) 

and {&m^) = {m\i^lt^\m) is the matrix element of the operator in the m-pion state. Excited states 
are suppressed in this expression and we assume ^ i^g"' ^ = for the cases we consider. For m- 
pion systems, contributions involving colour singlet sub-states propagating around the temporal 
boundary result in there being contributions from matrix elements of states with I <m pions. Near 
the middle of the temporal extent, t r^T /2, these contributions may be important. 



In the large T limit, only the £ = state contributes in Eqs. ( |2.2[ ) and ( [2.4[ ) and it is clear that 
the ratio 

= p^\..p^^{^'),nM, (2.5) 

where = {0,Em), This ratio will be independent of the sink and operator insertion times, t and 
z, provided fo ^ |f — t| <C T and determines the in-medium Mellin moment. It follows that the 
double ratio = Rm^ /R^^'^ determines the ratio of moments in medium to those in free space 
that we are interested in up to a simple kinematic factor. This double ratio is independent of the 
renormalisation scale, n, obviating the need for calculating the coefficients necessary to match 
lattice operators to operators in the MS scheme. If t — to is not much less than T, the above ratios 
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Label a [fm] xT niji [MeV] nijiL nijiT Measurements 



F 


0.09 


28^ X 96 


320 


4.1 


14.0 


432 


CI 


0.12 


20^ X 64 


290 


3.7 


11.7 


1450 


C2 


0.12 


20^ X 64 


350 


4.4 


14.2 


837 


C3 


0.12 


20^ X 64 


490 


6.2 


19.9 


1000 


L 


0.12 


20^ X 64 


350 


6.3 


14.2 


250 



Table 1: Details of measurements and ensembles used in this calculation. For each data set, the remaining 
columns correspond to the lattice spacing, lattice dimensions, valence pion mass and number of configura- 
tions X number of sources on each configuration. 



will be contaminated by thermal contributions. However, Eqs. (|2.4| ) and (|2.2| ) can still be used to 
extract the Mellin moments 

In terms of a lattice calculation, the two and three point functions we require are complicated 
by the many Wick contractions that arise in multi-pion systems. We proceed by defining partly 
contracted objects 

n(f;p) = 52'^«(x,f;xo,fo)'5d(x,f;xo,fo), 

X 

nW(?,T;p) = £^''P-X(y,T;xo,fo)757^''°£»^' •■•£»''" V55„(x,f;y,T)5;^(x,?;xo,fo), 

where 5^ is the quark propagator of flavour q and we have used 75-hermiticity, S{x,y) = y=,S'^ {y,x)y=,, 
to reverse the arguments of propagators from sink to source and contraction on spin and colour 
indices is assumed. These objects are can be viewed as time-dependent 12 x 12 matrices in spin 
and colour space. Using the techniques developed in Refs. [Q, |5|, ^ ^, the required contractions 
can now be built in terms of spin-colour traces of products of these matrices that can be efficiently 
computed. As the multi-hadron correlators decay very rapidly with Euclidean time, high precision 
arithmetic is needed in these calculations for which we use the QD library [||]. 

Numerical details: Our calculations are performed using gauge configurations generated by 
the MILC collaboration [|] using the rooted-staggered formulation of quarks and the asqtad gauge 
action. One level of HYP smearing [ |l0| ] is applied to these configurations to reduce short distance 
fluctuations. The ensembles used in this study are shown in Table [l], where we also report the 
number of configurations used and the number of source locations used on each configuration (the 



L ensemble has been used solely for checks of volume dependence). Domain-wall [11, 12] valence 



quark propagators have been calculated from APE smeared [13] sources at various locations by 



the NPLQCD and LHP collaborations []1J, [I5|]. These are then APE smeared on the sink time- 
slice with fixed momentum and used as the source for the sequential propagators connecting to the 
operator generated using the same action. The source-sink separation in the three point correlation 
functions is chosen at fsep/« = {t — ta)/<^ G {16, 20, 24, 28, 32} on the CI, C2 and C3 ensembles 
and at 4ep/fl: ^ {24, 32, 48} on the F ensemble. The operator insertion time is varied over the entire 
lattice. 

3. Thermal contamination 

A major issue in the present calculations is the contribution of thermal states to two- and three- 



4 



Medium Effects in Parton Distributions 



William Detmold 




point correlations functions. Given that the temporal extent of the lattice is fixed, one needs to go 
to early Euclidean times in order to ensure that the true ground state of the system is dominating 
the signal. However at early time, one must be concerned about (forward-going) excitations that 
are damped at later Euclidean times and a careful analysis is required. By using the full form of 
the expected correlator, Eq. ( |2.2D , and performing fits to all 12 correlators, the energies and the 
multiple overlap factors Zm can be determined. Since we successively fit C\, C2, ■ ■ ■ , Cn, only a 
single energy is determined by each fit and to account for correlations we perform these fits using 
a bootstrap procedure. The results of our fits to the correlators on the C 1 ensemble are shown in 
Fig-i 

The extracted bootstrap sets of parameters can then be used to reconstruct the correlator in the 
Umit of infinite temporal extent by removing all thermal contributions. In Fig. ^ the ratios of these 
zero temperature extrapolation of the correlators to the original correlators are shown as a function 
of Euclidean time for the m = 1, . . . , 12 pion systems (again for the CI ensemble). It is clear from 
this figure that for source-sink separations beyond t 20, it will prove difficult to extract the m- 
pion matrix element from the corresponding three-point function as the desired contribution will be 
suppressed relative to the thermal contributions. This is also the case for the C3 ensemble (where 
nijcT ~ 20) for larger numbers of pions as shown in the right panel of the figure. 



4. Quark momentum fraction 

In order to perform extractions of the matrix elements, we use the bootstrap list of the m-pion 
energies and overlap factors obtained from the two-point functions and input them into the expected 



spectral decomposition of the three-point function, Eq. ([2.4|). Using multiple different source-sink 
separations, we then fit the parameters {^m^) to the three point data under the bootstrap procedure. 
In Fig. ^, we show preliminary results for the the dependence of {x),nK on the density of the pion 
gas for the C3 ensemble. Mild dependence on the density of pions is observed but further work 
remains to better quantify systematic effects and to study the dependence on quark masses and the 
continuum and infinite volume limits. 
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Figure 2: The ratio of the zero temperature reconstruction of the correlators for m= 1 , . . . , 12 pion systems 
to the fitted (finite T) correlators. Data are from the CI (left) and C3 (right) ensembles. When this ratio 
deviates from unity, thermal effects are making significant contributions to the correlator. 

1.6 p ^ ^ ^ \ ^ ^ ^ \ ^ ^ ^ \ ^ ^ ^ \ ^ ^ ^ 1— 



1.4 - 



— 1.2 - 



0,8 U . . . \ . . . \ . . . \ . . . \ . . . U 

2 4 6 8 10 12 

m 

Figure 3: Extracted ratio of pion momentum fraction in an m-pion system to that in a single pion, {x),nnl (x) k 
for the C3 ensemble as a function of the number of pions in the system. 



5. Discussion 

It is clear from our preliminary studies presented here that the investigation of multi-hadron 
matrix elements is a challenging task in lattice QCD. While the current investigations focus on 
multi-pion systems, similar techniques in principle allow access to matrix elements in nuclei. With 
the recent observation of bound light nuclei in quenched QCD [ |T6| , [l7| ] and QCD at unphysical 
quark masses [11, |19|, H, further studies are warranted. However, the large statistics needed in 



lattice studies of few-body nuclei make this a daunting task for the future. On a more positive note, 
the thermal contamination that has hampered then current investigations is particularly vexing as 
Em-e + E£ < E,„ for these systems and in light nuclei, where the binding per nuclei increases, this 
may not be such a significant issue. 
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